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1. Introduction 

It is a pleasure to be able to speak at a meeting in San Diego in honour 
of Bill Helton, through whose early papers, especially [31], I first became 
interested in applications of operator theory to engineering. I shall discuss 
a problem of Heltonian character: a hard problem in pure analysis, with 
immediate applications in control engineering, which can be addressed by 
operator-theoretic methods. Furthermore, the main advances I shall describe 
are based on some highly original ideas of Jim Aglcr, so that San Diego is 
the ideal place for my talk. 

The /i-synthesis problem is an interpolation problem for analytic ma- 
trix functions, a generalization of the classical problems of Nevanlinna-Pick, 
Caratheodory- Fejer and Nehari. The symbol p, denotes a type of cost function 
that generalizes the operator and H°° norms, and the /i-synthesis problem is 
to construct an analytic matrix function F on the unit disc satisfying a finite 
number of interpolation conditions and such that fi(F(X)) < 1 for |A| < 1. 
The precise definition of /z is in Section [4] below, but for most of the paper we 
need only a familiar special case of /i - the spectral radius of a square matrix 
A, which we denote by r(A). 
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The purpose of this lecture is to present some cases of the /i-synthesis 
problem that are amenable to analysis. I shall summarize some results that 
are scattered through a number of papers, mainly by Jim Agler and me but 
also several others of my collaborators, without attempting to survey all the 
literature on the topic. I shall also say a little about recent results of some 
specialists in several complex variables which bear on the matter and may 
lead to progress on other instances of /i-synthesis. 

Although the cases to be described here are too special to have sig- 
nificant practical applications, they do throw some light on the /i-synthesis 
problem. More concretely, the results below could be used to provide test data 
for existing numerical methods and to illuminate the phenomenon (known to 
engineers) of the numerical instability of some /i-synthesis problems. 

We are interested in citeria for /i-synthesis problems to be solvable. Here 
is an example. We denote by D and T the open unit disc and the unit circle 
respectively in the complex plane C. 

Theorem 1.1. Let X\, X 2 eD be distinct points, let Wi, W 2 be nonscalar 2x2 
matrices of spectral radius less than 1 and let Sj = tr Wj , pj = det Wj for 
j = 1,2. The following three statements are equivalent: 
(1) there exists an analytic function F : D — > C 2x2 such that 

F(Ai) = Wi, F(X 2 ) = W 2 

and 

r(F(X)) < 1 for all A G D; 

(2) 

Ai — A2 
1 - A 2 Ai ; 
(3) 

1- -1 2 

(2 - usj)(2 - ujsj) - {2up t - Si)(2upj - sj) > 

for all w£T. 

The paper is organised as follows. Section [2] contains the definition of 
the spectral Nevanlinna-Pick problem, sketches the ideas that led to Theo- 
rem [TTT]- reduction to the complex geommetry of the symmetrized bidisc G, 
the associated "magic functions" $ w and the calculation of the Caratheodory 
distance on G - and fills in the final details of the proof of Theorem 11.11 us- 
ing the results of [IT] . It also discusses ill-conditioning and the possibility of 
generalization of Theorem 11.11 In Section [3] there is an analogous solvability 
criterion for a variant of the spectral Nevanlinna-Pick problem in which the 
two interpolation points coalesce (Theorem l3.1[) . In Section|H besides the def- 
inition of n and /i-synthesis, there is some motivation and history. Important 
work by H. Bercovici, C. Foia§ and A. Tannenbaum is briefly described, as is 
Bill Helton's alternative approach to robust stabilization problems. In Sec- 
tion [5] we consider an instance of /^-synthesis other than the spectral radius. 



(s 2 pi - S!p 2 )uj 2 + 2(p 2 - px)co + si - s 2 

max 



weT (si - s 2 pi)uj 2 - 2(1 - pip 2 )tu + s 2 - sip 2 



< 
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Here we can only obtain a solvability criterion in two very special circum- 
stances (Theorems 15. II and 15. 2p . The paper concludes with some speculations 
in Section |U 

We shall denote the closed unit disc in the complex plane by A. 



2. The spectral Nevanlinna-Pick problem 

A particularly appealing special case of the /i-synthesis problem is the spec- 
tral Nevanlinna-Pick problem: 

Problem SNP Given distinct points Ai,...,A„ G B and k x k matrices 
Wi, . . . , W n , construct an analytic k x k matrix function F onU) such that 

F(X j ) = W j for j = 1,..., n (2.1) 

and 

r{F(X)) < 1 for all A £ D. (2.2) 

When k = 1 this is just the classical Nevanlinna-Pick problem, and it 
is well known that a suitable F exists if and only if a certain n x n matrix 
formed from the Xj and Wj is positive (this is Pick's Theorem). We should 
very much like to have a similarly elegant solvability criterion for the case 
that k > 1, but strenuous efforts by numerous mathematicians over three 
decades have failed to find one. 

About 15 years ago Jim Agler and I devised a new approach to the 
problem in the case k = 2 based on operator theory and a dash of several 
complex variables ([5] to [H]). Since interpolation of the eigenvalues fails, 
how about interpolation of the coefficients of the characteristic polynomials 
of the Wj , or in other words of the elementary symmetric functions of the 
eigenvalues? This thought brought us to the study of the complex geometry 
of a certain set r C C 2 , defined below. By this route we were able to analyse 
quite fully the simplest then-unsolved case of the spectral Nevanlinna-Pick 
problem: the case n = k = 2. For the purpose of engineering application this 
is a modest achievement, but it nevertheless constituted progress. It had the 
merit of revealing some unsuspected intricacies of the problem, and may yet 
lead to further discoveries. 

2.1. The symmetrized bidisc V 

We introduce the notation 



r = {(z + w, zw) : z,w G A}, 
G = {(z + w, zw) : z,w £ D}. 



(2.3) 
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F and G are called the closed and open symmetrized bidiscs respectively. Their 
importance lies in their relation to the sets 

E d = {A e C 2x2 : r(A) < 1}, 

E° d = {A G C 2x2 : r(A) < 1}. 

E and its interior E° are sometines called "spectral unit balls", though the 
terminology is misleading since they are not remotely ball-like, being un- 
bounded and non-convex. Observe that, for a 2 x 2 matrix A, 

AeS'S the zeros of the polynomial A 2 — tr AX + det A lie in A 

tr A = z + w, dot A = zw for some z, w G A. 

We thus have the following simple assertion. 

Proposition 2.1. For any A G C 2x2 

A G E i/ and on/y i/ (tr A, det A) G T, 

A G E° i/ and on/?/ «/ (tr A, det A) G G. 

Consequently, if F : D — > E is analytic and satisfies the equations (|2.1j) 

above, where k = 2, then /i = (trF, det-F) is an analytic map from D to T 
satisfying the interpolation conditions 

h(Xj) = (tr Wj, det Wj) for j = 1, . . . , n. (2.4) 

Let us assume that none of the target matrices Wj is a scalar multiple of the 
identity. On this hypothesis it is simple to show the converse |16) by similarity 
transformation of the Wj to companion form. 

Proposition 2.2. Let Ai, . . . , A„ be distinct points in D and let W\, . . . , W n be 

nonscalar 2x2 matrices. There exists an analytic map F : D — 5- C 2x2 such 
that equations (|2.1[) and (|2.2[) hold if and only if there exists an analytic map 
h : D — > r that satisfies the conditions (|2.4[) . 

We have therefore (in the case k = 2) reduced the given analytic inter- 
polation problem for E-valued functions to one for T-valued functions (the 
assumption on the Wj is harmless, since any constraint for which Wj is scalar 
may be removed by the standard process of Schur reduction) . 

Why is it an advance to replace E by F? For one thing, of the two sets, 
the geometry of T is considerably the less rebarbative. E is an unbounded, 
non-smooth 4-complex-dimensional set with spikes shooting off to infinity in 
many directions. T is somewhat better: it is compact and only 2-complex- 
dimensional, though T too is non-convex and not smoothly bounded. But 
the true reason that T is amenable to analysis is that there is a 1-parameter 
family of linear fractional functions, analytic on G, that has special properties 
vis-a-vis T. For u in the unit circle T we define 

M**) - (2-5) 

2 — uis 

We use the variables s and p to suggest "sum" and "product" . The $ w de- 
termine G in the following sense. 
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Proposition 2.3. For every uj € T, $ w mapsG analytically intol}. Conversely, 
if (s,p) £ C 2 is such that |<£> w (s,p)| < 1 for all u> £ T, then (s,p) £ G. 

Both statements can be derived from the identity 



\2-z-w\ 2 - \2zw ~z-w\ 2 = 2(1 - |z| 2 )|l - w\ 2 + 2(1 - \w\ 2 )\l - z\ 2 . 
See [TO Theorem 2.1] for details. 



There is an analogous statement for T, but there are some subtleties. 
For one thing 3?^ is undefined at (2w,w 2 ) £ T when wel 

Proposition 2.4. For every to £ T. $ w maps T \ {(2cl>, ui 2 )} analytically into 
A. Conversely, if (s,p) £ C 2 is smc/i i/iai |$ w (rs, r 2 p)| < 1 for all u> £ T anrf 
< r < 1 i/ien (s,p) e T. 

In the second statement of the proposition the parameter r is needed: 
it does not suffice that ^^(s,^)! < 1 for all 10 £ T (in the case that p = 1 the 
last statement is true if and only if s 6 K, whereas for (s,p) 6 T, of course 



We found the functions ^ by applying Agler's theory of families of 
operator tuples [SJ [3] . We studied the family J- of commuting pairs of op- 
erators for which T is a spectral set, and its dual cone J- 1 - (that is, the 
collection of hereditary polynomials that are positive on F) . Agler had pre- 
viously done the analogous analysis for the bidisc, and shown that the dual 
cone was generated by just two hereditary polynomials; this led to his cele- 
brated realization theorem for bounded analytic functions on the bidisc. On 
incorporating symmetry into the analysis we found that the cone J- 1 - had 
the 1-parameter family of generators 1 — <&^$ w , w£T. From this fact many 
conclusions follow: see [13] for more on these ideas. 

Operator theory played an essential role in our discovery of the functions 
Once they are known, however, the geometry of G and T can be developed 
without the use of operator theory. 

2.2. A necessary condition 

Suppose that F is a solution of the spectral Nevanlinna-Pick problem (|2.ip . 
(|2.2p with k = 2. Let us write Sj — tr Wj, Pj — det Wj for j = 1, . . . , n. For 
any w£T and < t < 1 the composition 



1*1 < 2). 




is an analytic self-map of D under which 



, , 9 n 2wi 2 p 

<\j !-> &u> {tSj , t Pj ) = 




for j = 1, . . . , 



n. 



Thus, by Pick's Theorem, 



1 - <S> ul {ts il t 2 p i )<S> ul {ts 0l t 2 p ] ) 
1 - AjAj 



> 0. 



(2.6) 
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E satisfying the 



On conjugating this matrix inequality by diag{2 — ujtsj} and letting a = tuj 
we obtain the following necessary condition for the solvability of a 2 x 2 
spectral Nevanlinna-Pick condition [5J Theorem 5.2]. 

Theorem 2.5. If there exists an analytic map F : B 
equations 

F(X j ) = W j for j ! n 

and 

r(F(A)) < 1 for allXeB 
then, for every a such that \a\ < 1, 



whe 



(2 - as l ){2 - asj) - \a\ z {2ap l - s i )(2ap j - sj) 
1 — AiA,' 



tr Wj , pj = det Wj for j = 1, . 



> 



(2.7) 



In the case that the Wj all have spectral radius strictly less than one, 
the condition (|2.7p holds for all a £ A if and only if it holds for all a £ T, and 
hence the condition only needs to be checked for a one-parameter pencil of 
matrices. It is of course less simple than the classical Pick condition in that 
it comprises an infinite collection of algebraic inequalities, but it is neverthe- 
less checkable in practice with the aid of standard numerical packages. Its 
major drawback is that it is not sufficient for solvability of the 2x2 spectral 
Nevanlinna-Pick problem. 

Example 2.6. Let < r < 1 and let 

A 2 A(A 3 +r)' 



^n 2(1 - r) l + rA3< 1 + .A3 

Let Ai,A2,A3 be any three distinct points in D and let h(\j) = (sj,Pj) for 
j = 1, 2, 3. We can prove [5] that, in any neighbourhood of (si, S2, S3) in (2B) 3 , 
there exists a point (s'i,s' 2 ,s' 3 ) such that (s'j,pj) G G, the Nevanlinna-Pick 
data 

A. • > -I'.. ; ;^./',K j = 1,2,3, 
are solvable for all ui £ T, but the Nevanlinna-Pick data 

A 3 h-> $ m ( Aj )(s^-Pi)> 3 = 1,2,3, 
are unsolvable for some Blaschke factor m. It follows that the interpolation 
data 

Xj^(s'j,Pj), j = 1,2,3, 
satisfy the necessary condition of Theorem 12.51 for solvability, and yet there 



is no analytic function h : B — >• T such that h(Xj) = (s'j,pj) for j = 1, 2, 3. 

Hence, if we choose nonscalar 2x2 matrices Wi,W2,Wa such that 
(tr Wj, det Wj) — (sj,pj), then the spectral Nevanlinna-Pick problem with 
data Xj i-> Wj satisfies the necessary condition of Theorem 12.51 and yet has 
no solution. 



See also [32] for another example. 
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2.3. Two points and two-by-two matrices 

When n = k = 2 the condition in Theorem [33J is sufficient for the solvability 
of the spectral Nevanlinna-Pick problem. 

We shall now prove the main theorem from Section [T] Recall the state- 
ment: 

Theorem 1.1. Let Ai, A2 S D be distinct points, let Wt, W 2 be nonscalar 2x2 
matrices of spectral radius less than 1 and let Sj = tr Wj , Pj = det Wj for 
j = 1,2. The following three statements are equivalent: 

(1) there exists an analytic function F : D — > C 2x2 such that 
F(Ax) = Wi, F{\ 2 ) = W 2 



and 



(2) 



max 



r(^(A)) < 1 for all A e ID 
(s 2 pi - S1P2V 2 + 2(p 2 - Pi)w + si - s 2 



(si - s 2 pi)u} 2 - 2(1 - pip 2 )u + s 2 - sip 2 



< 



Ai — A2 



(3) 



(2 — u)Si)(2 — uJSj) — (2ujpi — Si)(2upj — Sj) 

1 — A; A, 



1 - A2A1 



> 



(2.8) 



(2.9) 



for all w£T. 

The proof depends on some elementary notions from the theory of invariant 
distances. A good source for the general theory is [35], but here we only need 
the following rudiments. 

We denote by d the pseudohyperbolic distance on the unit disc D: 

Ai — A2 



d(Ai,A a ) = 



for Ai, A 2 £ 



I-A2A1 

For any domain G C" we define the Lempert function 5q : 51 x Q — > K + by 

5n(z u z 2 ) =infd(Ai,A 3 ) (2.10) 

over all Ai,A2 eD such that there exists an analytic map h : ED — > ft such 
that h(X\) — Z\ and h(X 2 ) = z 2 . We defin^l the Caratheodory distance Cn : 
ft x O -> K+ by 

C n (zi,z 2 ) =supd(/(2;i),/(«2)) (2.11) 
over all analytic maps / : — > B. If f2 is bounded then Cq is a metric on 57. 

It is not hard to see (by the Schwarz-Pick Lemma) that Co < 8q for 
any domain fl. The two quantities Co, Sq are not always equal - the punc- 
tured disc provides an example of inequality. The question of determining 
the domains Jl for which Co = ^o is one of the concerns of invariant distance 
theory. 



1 Conventionally the definition of the Caratheodory distance contains a tanh 1 on the right 
hand side of 1 12. lit . For present purposes it is convenient to omit the tanh -1 . 
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Proof. Let zj — (sj,Pj) £ G. 

(l)<=>-(2) In view of Proposition 12.21 we must show that the inequality 
is equivalent to the existence of an analytic h : D — >• T such that h(Xj) = zj 
for j — 1,2. By definition of the Lempert function 6&, such an h exists if and 
only if 

5 G (zi,z 2 ) < d(zi,z 2 ). 

By jTTJ Corollary 5.7] we have 5& = Cq, and by [TTJ Theorem 1.1 and Corol- 
lary 3.4], 



C G (z 1 ,z 2 ) 



maxd($ u (zi),$ u (z 2 )) 



max 



(s 2 pi - sxp 2 )uj 2 + 2(p 2 - pi)ui + si 



s 2 



(si - s 2 pi)uj 2 - 2(1 - pip 2 )uj + s 2 - Sip 2 



Thus the desired function h exists if and only if the inequality 



(2.12) 



holds. 



(2)<^>(3) By equation (|2.12|) . the inequality 
d(<Z>Ml),$M2))<d(\i,\ 2 



is equivalent to 
for all uj ET. 



By the Schwarz-Pick Lemma, this inequality holds if and only if, for all uj G T, 
there exists a function f u in the Schur class such that / w (A_y) = &w(zj) for 
j = 1, 2. By Pick's Theorem this in turn is equivalent to the relation 

, ~i 2 

> 0. 



$ w (Zi)$u,(*j) 



1 



Conjugate by diag{2 — ujsi, 2 — WS2} to obtain (2)<=>(3). 



□ 



Remark 2.7. If one removes the hypothesis that W\,W 2 be nonscalar from 
Theorem 11.11 one can still give a solvability criterion. If both of the Wj 
are scalar matrices then the problem reduces to a scalar Nevanlinna-Pick 
problem. If Wi = cl and W 2 is nonscalar then the corresponding spectral 
Nevanlinna-Pick problem is solvable if and only if 

r{{W 2 - d){I - cW 2 )- 1 ) <d(Ai,A 2 ) 

(see [7J Theorem 2.4]). This inequality can also be expressed as a somewhat 
cumbersome algebraic inequality in c, s 2 ,p 2 and d(X\, X 2 ) Theorem 2.5(2)]. 

2.4. Ill-conditioned problems 

The results of the preceding subsection suggest that solvability of spectral 
Nevanlinna-Pick problems depends on the derogatory structure of the target 
matrices - that is, in the case of 2 x 2 matrices, on whether or not they 
are scalar matrices. It is indeed so, and in consequence problems in which a 
target matrix is close to scalar can be very ill-conditioned. 



Example 2.8. 



[3 Example 2.3] Let (3 e 
W l (a)= 0Q 



\ {0} and, for a e C let 

p 
ni 

U 1+/3 
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Consider the spectral Nevanlinna-Pick problem with data i-> Wi(a), f3 H> 
W-2,. If a — then the problem is not solvable. If a ^ 0, however, by Proposi- 
tion the problem is solvable if and only if there exists an analytic function 
/ : D -)■ T such that 

/(0) = (0,0) and/(/?) = T ^. 
It may be checked [5] that 

'2(1-0)* A(A-/3) 



/(A) V ' 1-/3A 

is such a function. Thus the problem has a solution F a for any a^O. Consider 
a sequence (a„) of nonzero complex numbers tending to zero: the functions 
F an cannot be locally bounded, else they would have a cluster point, which 
would solve the problem for a = 0. If a is, say, 10~ 100 then any numeri- 
cal method for the spectral Nevanlinna-Pick problem is liable to run into 
difficulty in this example. 

2.5. Uniqueness and the construction of interpolating functions 

Problem SNP never has a unique solution. If F is a solution of Problem 
SNP then so is p- 1 FP for any analytic function P : D -> C fexfe such that 
P(A) is nonsingular for every A £ D and P(\j) is a scalar matrix for each 
interpolation point Xj. There are always many such P that do not commute 
with F, save in the trivial case that F is scalar. Nevertheless, the solution of 
the corresponding interpolation problem for T can be unique. Consider again 
the case n = k = 2 with W\,Wi nonscalar. By Theorem ll.lt the problem 
is solvable if and only if inequality (|2.8| holds. In fact it is solvable uniquely 
if and only if inequality (|2.8[) holds with equality. This amounts to saying 
that each pair of distinct points of G lies on a unique complex geodesic of 
G, which is true by [13 Theorem 0.3]. (An analytic function h : D — > G is 
a complex geodesic of G if h has an analytic left-inverse). Moreover, in this 
case the unique analytic function h : D — > G such that h(Xj) = (sj,Pj) for 
j = 1,2 can be calculated explicitly as follows [HI Theorem 5.6]. 

Choose an wo G T such that the maximum on the left hand side of (|2.8j) 
is attained at ujq. Since equality holds in (|2.8j) . we have 



d($ Uo (zi),$ Uo (z 2 )) = d(Ai,A 2 ), 

where Zj — (sj,Pj). Thus <& Uo is a Caratheodory extremal function for the 
pair of points Z\, Z2 in G. It is easy (for example, by Schur reduction) to find 
the unique Blaschke product p of degree at most 2 such that 

p(Ai)=Pi, p(A 2 )=p 2 and p(u )=wl). 

Define s by 

s(A) = 2^ A )~ A for A E D. 

1 — UJqX 

Then h ^ (s,p) is the required complex geodesic. 
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Note that h is a rational function of degree at most 2. It can also be 
expressed in the form of a realization: h(X) — (tr H(X), det H(\)) where H is 
a 2 x 2 function in the Schur class given by 

H{\) = D + CX{\- A\y 1 B 



for a suitable unitary 3 x 3 or 4 x 4 matrix 
(see 0], [12 Theorem 1.7]). 



AB 
CD 



given by explicit formulae 



2.6. More points and bigger matrices 

Our hope in addressing the case n = k = 2 of the spectral Nevanlinna-Pick 
problem was of course that we could progress to the general case. Alas, we 
have not so far managed to do so. We have some hope of giving a good 
solvability criterion for the case k = 2, n = 3, but even the case n = 4 
appears to be too complicated for our present methods. 

The case of two points and k x k matrices, for any k, looks at first 
sight more promising. There is an obvious way to generalize the symmetrized 
bidisc: we define the open symmetrized polydisc Gfc to be the domain 

G k = {( CTl (z), . . . ,a k (z)) : z 6 B fe } C C k 

where a m denotes the elementary symmetric polynomial in z = (z 1 , . . . , z ) 
for 1 < m < k. Similarly one defines the closed symmetrized polydisc Tfe. As 
in the case k = 2, one can reduce Problem SNP to an interpolation problem 
for functions from D to under mild hypotheses on the target matrices Wj 
(specifically, that they be nonderogatory) . However, the connection between 
Problem SNP and the corresponding interpolation problems for T k are more 
complicated for k > 2, because there are more possibilities for the rational 
canonical forms of the target matrices [37] • The analogues for F k of the <E> W 
were described by D. J. Ogle [35] and subsequently other authors, e.g. [53J 
[2"5j . Ogle generalized to higher dimensions the operator-theoretic method of 
[B] and thereby obtained a necessary condition for solvability analogous to 
Theorem HH 

The solvability of Problem SNP when n = 2 is generically equivalent to 
the inequality 

SG k {zi,z 2 ) < d(Xi,X 2 ) 
where Zj is the fctuple of coefficients in the characteristic polynomial of Wj . 
All we need is an effective formula for <5c fc ■ It turns out that this is a much 
harder problem for k > 2. In particular, it is false that 5g, k = C*c fc when 
k > 2. This discovery [35] was disappointing, but not altogether surprising. 

There is another type of solvability criterion for the 2x2 spectral 
Nevanlinna-Pick problem with general n [TU1 Q3], but it involves a search 
over a nonconvex set, and so does not count for the purpose of this paper as 
an analytic solution of the problem. Another paper on the topic is [24) . 

It is heartening that the study of the complex geometry and analysis 
of the symmetrized polydisc has been taken up by a number of specialists 
in several complex variables, including G. Bharali, C. Costara, A. Edigarian, 
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M. Jarnicki, L. Kosinski, N. Nikolov, P. Pflug, P. Thomas and W. Zwonek. 
Between them they have made many interesting discoveries about these and 
related domains. There is every hope that some of their results will throw 
further light on the spectral Nevanlinna-Pick problem. 



3. The spectral Caratheodory-Fejer problem 

This is the problem that arises from the spectral Nevanlinna-Pick problem 
when the interpolation points coalesce at 0. 

Problem SCF Given kxk matrices Vq, V\, . . . , V n , find an analytic function 
F : D -> C kxk such that 

FU)(0) = V j for j = 0,..., n (3.1) 

and 

r(F(X)) < 1 for all A G D. (3.2) 



This problem also can be converted to an interpolation problem for ana- 
lytic functions from D into [HH Theorem 2.1], [37]. However, the resulting 
problem is again hard when k > 2, and the only truly explicit solution we 
have is in the case k = 2, n = 1 34, Theorem 1.1]. 



Theorem 3.1. Let 



V m = k" 112 



for m = 0, 1 



and suppose that Vq is nonscalar. There exists an analytic function F 
C 2x2 such that 



F(0) = V , 
if and only if 



F'(0) = Vx and r(F(X)) < 1 for all A G 



max 

|U)|=1 



(sxpo - s pi)uj 2 + 2ujpi 



91 



u) 2 (s ~ s po) - 2w(l - |po| 2 ) + s - s po 



< 1, 



(3.3) 



(3.4) 



where 



s =trVb, po = detVb, 
si = trVi, pi 



"11 V 12 
v 21 v 22 



Vn-i Wot 



The proof of this theorem in 34 again depends on the calculation in 
|llj of the Caratheodory metric on G, but this time on the infinitesimal 
version cq of the metric: the left hand side of inequality (|3.4p is the value 
of cg at (so,Po) m the direction (si,pi). This fact is [HJ Corollary 4.4], but 
unfortunately there is an oj missing in the statement of Corollary 4.4. The 
proof shows that the correct formula is as in ()3.4j) . An important step is 
the proof that the infinitesimal Caratheodory and Kobayashi metrics on G 
coincide. 

The ideas behind Theorem 13 . 1 1 can be used to find solutions of Problem 
SCF: see [MJ Section 6]. The ideas can also be used to derive a necessary 
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condition for the spectral Caratheodory-Fejer problem (|3.1[) . (13.21) in the case 
that n = 1 and fc > 2 [3U Theorem 4.1], but there is no reason to expect this 
condition to be sufficient. 



4. The structured singular value 

The structured singular value of a matrix relative to a space of matrices was 
introduced by J. C. Doyle and G. Stein in the early 1980s [25] and was 
denoted by /i. It is a refinement of the usual operator norm of a matrix 
and is motivated by the problem of the robust stabilization of a plant that is 
subject to structured uncertainty. Initially, in the H°° approach to robustness, 
the uncertainty of a plant was modelled by a meromorphic matrix function 
(on a disc or half plane) that is subject to an L°° bound but is otherwise 
completely unknown. The problem of the simultaneous stabilization of the 
resulting collection of plant models could then be reduced to some classical 
analysis and operator theory, notably to the far-reaching results of Adamyan, 
Arov and Krein from the 1970s 30 . 

In practice one may have some structural information about the uncer- 
tainty in a plant - for example, that certain entries are zero. By incorporating 
such structural information one should be able to achieve a less conservative 
stabilizing controller. The structured singular value was devised for this pur- 
pose. A good account of these notions is in 27, Chapter 8]. Unfortunately, 
the behaviour of /i differs radically from that of the operator norm - for one 
thing, /i is not in general a norm at all, and none of the relevant classical 
theorems (such as Pick's theorem) or methods appear to extend to the cor- 
responding questions for \x. This provides a challenge for mathematicians: 
we should help out our colleagues in engineering by creating an AAK-type 
theory for /i. 

For any A S C kxl and any subspace E of C txh we define the structured 
singular value ^e(A) by 

— ^— = inf{||A|| : X 6 E, 1 - AX is singular} (4.1) 
Ve{A) 

with the understanding that /j,e(A) = if 1 — AX is always nonsingular. 

Two instances of the structured singular value are the operator norm 
||.|| (relative to the Euclidean norms on C fc and C £ ) and the spectral radius 
r. If we take E = <C txk then we find that (1e{A) = \\A\\. On the other hand, 
if k = I and we choose E to be the space of scalar multiples of the identity 
matrix, then /j,e(A) = r(A). These two special [is are in a sense extremal: it 
is always the case, for any E, that he{A) < \\A\\. If k = £ and E contains 
the identity matrix, then [Xe{A) > r(A). A comprehensive discussion of the 
properties of /x can be found in [40 . 

Here is a formulation of the /j,- synthesis problem |26[ I27j . 

Given positive integers k, I, a subspace E ofC ixk and analytic functions 
A, B,C on D of types k x £,k x k and £ x I respectively, construct an analytic 
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function F : D — > C kxe of the form 

F = A + BQC for some analytic Q:B^C kxl (4.2) 

such that 

Ve{F(X)) < 1 for all A G D. (4.3) 

The condition (|4.2|) . that F be expressible in the form A + BQC for 
some analytic Q, can be regarded as an interpolation condition on F. In the 
event that k = £ , B is the scalar polynomial 

B(X) = (X-X 1 )...(X-X n )I 

with distinct zeros A^ £ B and C is constant and equal to the identity, then 
F is expressible in the form A + BQC if and only if 

F(X 1 )=A(X 1 ),...,F(X n )=A(X n ). 

With this choice of B and C, if we take E to be the space of scalar matrices, 
we obtain precisely the spectral Nevanlinna-Pick problem. If we now replace 
B by the polynomial A™, we get the spectral Caratheodory-Fejer problem. 

In engineering applications /i-synthcsis problems arise after some anal- 
ysis is carried out on the plant model to produce the A, B and C in condition 
(|4.2|) . and the resulting B and C will not usually be scalar functions. Never- 
theless, explicit pointwise interpolation conditions provide a class of easily- 
formulated test cases, and it is arguable that such problems are the hardest 
cases of //-synthesis. 

Conditions of the form (|4.2p are said to be of model matching type [30] . 

The most sustained attempt to develop an AAK-type theory for the 
structured singular value in full generality is due to H. Bercovici, C. Foias, and 
A. Tannenbaum f |15] to|21]V They have a far-reaching theory: inter alia they 
have constructed many illuminating examples, found properties of extremal 
solutions and obtained a type of solvability criterion for /i-synthesis problems. 
The criterion results from a combination of the Commutant Lifting Theorem 
with the application of similarity transformations. To apply the criterion to 
a concrete spectral Nevanlinna-Pick problem one must solve an optimization 
problem over a high-dimensional unbounded and non-convex set. We can 
certainly hope that this is not the last word on the subject of solvability. 
Despite the achievements of Bercovici, Foia§ and Tannenbaum, there is still 
plenty of room for further study of /z-synthesis. 

One of their examples [18| Section 7, Example 5] exhibits an important 
fact about the spectral Nevanlinna-Pick problem: diagonalization does not 
work. It shows that diagonalization of the target matrices Wj in Problem 
SNP by similarity transformations, even when possible, does not help solve 
the problem. One could hope that if the Wj were diagonal one might be able 
to decouple the problem into a series of scalar interpolation problems, but 
they show that such a hope is vain. 

Bill Helton himself, along with collaborators, has developed an alter- 
native approach to the refinement of H°° control; his viewpoint is set out 
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in |32] . His part in the introduction of the results of Adamyan, Arov, Krein 
and other operator-theorists into robust control theory in the early 1980s 
is well known. He subsequently worked extensively (with Orlando Merino, 
Trent Walker and others) during the 1990s on the more delicate optimization 
problems that arise from refinements of the basic H°° picture of modelling 
uncertainty. As in the fi approach, the aim is to incorporate more subtle 
specifications and robustness conditions into methods for controller design. 
He developed a very flexible formulation of such problems as optimization 
problems over spaces of vector- valued analytic functions on the disc, and de- 
vised an algorithm for their numerical solution - see |33j and several other 
papers. The authors proved convergence results and described numerical tri- 
als. However, the spectral Nevanlinna-Pick problem cannot be satisfactorily 
treated by the Helton scheme. Although it can be cast in the basic problem 
formulation [32 Chapter 2], solution algorithms require smoothness proper- 
ties (of the function £ T") which the spectral radius does not possess. 



5. The next case of /i 

After the two extremes fi — and fi — r the next natural case to 

consider is the one in which, in ()4.1|) . k — £ and E is the space Diag(fc) of 
diagonal matrices. For the rest of this section fi will denote /^Diag(2) and we 
shall study the following problem: 

Given distinct points Ai, . . . , A n £ D and 2x2 matrices Wi, . . . , W n , 
construct an analytic function F : D — > C 2x2 such that 

F(X j ) = W j forj = l,...,n (5.1) 

and 

fj,(F(\)) < 1 for all A e D. (5.2) 



For the 2x2 spectral Nevanlinna-Pick problem we had some modest 
success through reduction to an interpolation problem for T-valued functions. 
In the present case we tried an analogous approach, with still more modest 
success [TJ [2j [4TJ . The following result is [21 Theorem 9.4 and Remark 9.5(iii)]. 

Theorem 5.1. Let A £ D, A ^ 0, let ( e C and let 



Wi 



0C 




w 2 



a * 
* b 



(5.3) 



Suppose that \b\ < \a\ and let p — detW^- There exists an analytic function 
F : B — > C 2x2 such that 



F{X 1 ) = W 1 , F(\ 2 ) = W 2 and /z(F(A)) < 1 for all A £ D (5.4) 
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(|a| 2 + \b\ 2 + 2\ab -p|)|A | 2 + \p\ 2 > if C = 0. 



The stars in the formula for W2 in (|5.3p denote arbitrary complex num- 
bers. 

What is the analog of T for this case of /z? To determine whether a 2 x 2 
matrix A — [ay] satisfies r(A) < 1 one needs to know only the two numbers 
tr A and det A; this fact means that the spectral Nevanlinna-Pick problem 
can generically be reduced to an interpolation problem for T. To determine 
whether fi(A) < 1 one needs to know the three numbers an, 022, det A. This 
led us to introduce a domain E which we call the tetrablock: 

E = {x eC 3 : 1 - x 1 z - x 2 w + x 3 zw ^ whenever \z\ < 1, \w\ < 1}. (5.5) 

Its closure is denoted by E. The name reflects the fact that the intersection 
of E with M 3 is a regular tetrahedron. The domain E is relevant because 
n(A) < 1 if and only if (an, 022, det A) G E. There exists a solution of the 2- 
point /i-synthesis problem (I5.4[) if and only if the corresponding interpolation 
problem for analytic functions from D to E is solvable O Theorem 9.2], and 
accordingly the solvability problem for this //-synthesis problem is equivalent 
to the calculation of the Lempert function Sr. As far as I know no one has yet 
computed Se for a general pair of points of E, but we did calculate it in the 
case that one of the points is the origin in C 3 , that is, we proved a Schwarz 
lemma for E. The result is Theorem 15. II 

Observe that ill-conditioning appears in this instance of /i-synthesis too 
Remark 9.5(iv)]. If, in Theorem 15. 1[ a = b = p = i then there exists a 
solution Fq of the problem if and only if 



sensitive to small changes in £ near 0. 

The complex geometry of E has also proved to be of interest to re- 
searchers in several complex variables. To my surprise, it was recently shown 
[2"S] that the Lempert function and the Caratheodory distance on E coincide. 
This might be a step on the way to the derivation of a formula for i5e- It 
would suffice to compute Se in the case that one of the two points is of the 
form (0, 0, A) for some A € [0, 1), since every point of E is the image of such 
a point under an automorphism of E pfl"} Theorem 5.2]. 

The fourth and final special case of //-synthesis in this paper is the 
//-analog of the 2x2 Caratheodory- Fejer problem: 




Thus if I < |Aq| < -4=, the are not locally bounded as C —> 0, and so are 
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Given 2x2 matrices Vq, . . . , V n , construct an analytic function F : B — > 
C 2x2 such that 

F (j) (0) = Vj forj = 0,...,n and /z(F(A)) < 1 for all AeD. 

Again the problem can be reduced to an interpolation problem for E, 
but the resulting problem has only been solved in an exceedingly special case. 

Theorem 5.2. Let Vo, V\ be 2 x 2 matrices such that 

Vq ~ [oo 

for some ( £ C and V% — [vij] is nondiagonal. There exists an analytic 
function F : D — > C 2x2 such that 

F(0) = V , F'(0) = Vi and fi{F(X)) < 1 for all AeD 

if and only if 

max{|wn|, 1^22 1} + |C«2i| < 1- 
This result follows from [HI Theorem 2.1]. 



6. Conclusion 

Although /i-analysis remains a useful tool, it is fair to say that /i-synthesis, 
as a major technique for robust control system design, has been something of 
a disappointment up to now. The trouble is that the /z-synthesis problem is 
difficult. It is a highly non-convex problem. There do exist heuristic numerical 
methods for addressing particular /z-synthesis problems, notably a Matlab 
toolbox [33] based on the U DK algorithm" [27J Section 9.3], but there is no 
practical solvability criterion, no fast algorithm nor any convergence theorem 
for any known algorithm. For these reasons engineers have largely turned to 
other approaches to robust stabilization over the past 20 years. If, however, a 
satisfactory analytic theory of the problem is developed, engineers' attention 
may well return to zz-synthesis as a promising design tool. We are still far 
from having such a theory, but perhaps these special cases and the interest 
of the several complex variables community may yet lead to one. 
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